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Abstract 

Much  of  seismic  nuclear  monitoring  depends  on  a  proper  understanding  of  attenuation 
and  scattering  of  regional  seismic  waves.  This  report  seeks  to  characterize  wave  scattering 
from  both  an  observational  and  theoretical  approach.  A  theory  of  wave  propagation  in 
fractally  homogeneous  random  media  is  developed  and  applied  to  southern  Tibet  in  order 
to  quantify  the  analytical  and  statistical  properties  of  crust  and  upper  lithosphere  hetero¬ 
geneity.  Seismic  velocity  heterogeneity  is  modeled  by  a  self-affine  band  limited  random 
fractal.  An  approximation  to  multiple  scattering  theory  yields  a  relationship  between  the 
exponent  of  the  frequency  dependent  scattering  quality  factor  Q  and  the  fractal  dimension 
of  velocity  iso-surfaces.  In  this  study,  the  attenuation  quality  factor  of  P-waves,  Qp,  is 
calculated  for  southern  Tibet  and  the  implications  for  the  fractal  dimension  of  the  veloc¬ 
ity  field  are  discussed.  In  addition,  the  topic  of  intermittent  heterogeneity  is  introduced. 
The  generalization  of  the  model  to  fractally  homogeneous  media  implies  that  scattering 
attenuation  measurements  are  characterized  by  two  fractal  dimensions:  one  related  to  the 
self-affinity  of  the  velocity  field  and  the  other  related  to  the  metric  of  the  space  which  the 
field  occupies. 

Within  a  fractal  random  media  model,  the  attenuation  calculation  for  southern  Tibet 
yields  a  fractal  dimension  of  velocity  iso-surfaces  of  =  9/4  which  implies  the  existence 
of  complex  large-scale  structure  in  the  upper  lithosphere.  For  the  more  general  case  of 
fractally  homogeneous  heterogeneity,  the  frequency  dependence  of  the  scattering  quality 
factor  is  a  function  of  two  parameters  D  and  another  fractal  dimension  D.  The  predicted 
frequency  exponent  can  differ  from  the  homogeneous  case  by  as  much  as  2/3  for  reason¬ 
able  choices  of  the  parameters  D  and  D.  However,  the  relative  contributions  of  topological 
disorder  and  intermittency  are  inextricably  linked  and  cannot  be  isolated  by  scattering 
attenuation  measurements.  Continued  investigation  of  wave  propagation  in  fractally  ho¬ 
mogeneous  media  may  further  quantify  the  statistical  properties  of  crustal  heterogeneity. 
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Objectives 

Much  of  seismic  nuclear  monitoring  depends  on  an  understanding  of  the  attenuation 
and  scattering  characteristics  of  regional  seismic  phases.  The  objective  of  this  report  is 
to  quantify  the  nature  of  heterogeneity  in  the  crust  and  upper  lithosphere  in  southern 
Tibet  and  to  show  that  for  wave  propagation  in  general  fractally  homogeneous  media,  the 
frequency  dependence  of  the  scattering  attenuation  quality  factor  Q  is  a  function  of  two 
parameters  related  to  the  geometry  of  the  velocity  field.  Following  Wu  and  Aki,  seismic 
velocity  heterogeneity  is  modeled  by  a  band  limited  random  fractal  (Wu  and  Aki,  1985).  An 
approximation  to  multiple  scattering  theory  in  such  media  yields  a  relationship  between  the 
exponent  of  the  frequency  dependent  quality  factor  and  the  fractal  dimension  of  velocity 
iso-surfaces  (Wu,  1982).  In  this  study,  the  attenuation  quality  factor  of  P-waves,  Qp,  is 
calculated  for  southern  Tibet  and  the  implications  for  the  fractal  dimension  of  the  velocity 
field  are  discussed.  In  addition,  the  topic  of  intermittent  heterogeneity  is  introduced. 
The  generalization  of  the  model  to  fractally  homogeneous  media  implies  that  scattering 
attenuation  measurements  are  characterized  by  two  fractal  dimensions:  one  related  to  the 
self-afl&nity  of  the  velocity  field  and  the  other  related  to  the  metric  of  the  space  which  the 
field  occupies. 


Research  Accomplished 

Model  Random  Media 

The  model  random  media  is  defined  by  a  position  dependent  velocity  c{f)  =  cq  -|- 
6c(r)  where  the  fluctuating  part  Sc{f}  is  a  zero-mean,  isotropic,  random  funtion.  For 
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the  purposes  of  this  study  the  fluctuation  is  taken  to  be  a  fractional  Brownian  function 
8c{f)  =  (Mandelbrot,  1977).  The  fractional  Brownian  function  Bnir)  is  defined 

through  a  variance  or  structure  function, 

F(r)  =<  {Bnix'  +  (1) 


For  the  value  of  the  parameter  H  —  1/2,  the  function  Bh{t)  along  any  line  in  space  is 
a  Brownian  walk  whose  variance  scales  linearly  with  sample  size.  Mandelbrot  has  shown 
that  the  fractal  dimension  of  the  iso-surfaces  of  such  a  field  is  D  =  S  —  H  (Mandelbrot, 
1975).  Thus,  if  the  fluctuations  are  characterized  by  a  root  mean  square  variation  7,  the 
structure  function  can  be  written  as  F{r)  —  The  value  of  F{r)  characterizes 

the  intensity  of  fluctuations  with  length  scale  smaller  than  or  on  the  order  of  r.  The 
one- dimensional  spectral  density  is  related  to  the  structure  function  through  an  integral 
transform  (Tatarski,  1960).  One  finds  that  the  spectral  density  is  given  by. 


$(fc)  = 


r(l  -h  2(3  -  D))  sin(7r(3  -  D))  ^  (2(3-d)+i) 

27r  ^  ^ 


(2) 


Of  course  any  real  medium  posseses  some  outer  scale  of  heterogeneity  such  that  the  power 
law  spectrum  is  truncated  at  some  point.  A  spectrum  that  corresponds  to  Eq.  (2)  in  the 
high  frequency  limit  and  includes  a  characteristic  length  is  the  von  Karman  spectrum, 


ilk)  =  7^r(3  -D  +  1/2) 

'  ’  (l  +  it2a2)(D-3+i/2),(3) 


where  a  is  the  outer  length  scale  of  the  inhomogeneity.  Wu  (Wu,  1982)  has  shown  that 
within  a  cumulative-forward,  single-backscatter  approximation  to  multiple  scattering  the¬ 
ory,  the  inverse  of  the  scattering  quality  factor 

g-'  =  2komV2ko)  -  $(2A:o)},  (4) 


where  ko  is  the  wavenumber  in  a  homogeneous  medium  of  velocity  Cq.  Thus,  substituting 
for  the  spectral  density  $,  in  the  high  frequency  regime  the  model  yields  a  frequency 
dependent  quality  factor, 

g-i  (5) 

If  the  fractal  dimension  of  constant  velocity  surfaces  D  lies  between  2  <  D  <  3,  then 
Q~^  ~  k~^  where  0  <  <  2.  Thus,  observations  of  smaller  ^  correspond  to  larger  ampli¬ 

tudes  of  high  frequency  spectral  harmonics  and  smaller  amplitudes  of  low  frequency  com¬ 
ponents.  This  implies  the  existence  of  small  scale  field  structure  and  the  absence  of  large 
scale  structure.  This  is  illustrated  in  Figure  1  for  two  approximations  to  two-dimensional 
fractional  Brownian  functions.  As  the  value  of  D  increases  from  two  to  three  the  surfaces 
become  more  erratic  and  small  scale  detail  is  increased.  This  fact  is  the  underlying  geo¬ 
metrical  basis  for  the  observed  correlation  between  degree  of  frequency  dependence  of  Q 
and  crustal  stability  (Aki,  1980). 
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Southern  Tibet 


A  three-component  broadband  array  in  southern  Tibet  was  deployed  from  May  2,  1994 
to  October  23,  1994.  The  INDEPTH- II  passive  source  experiment  used  both  short  period 
and  broadband  seismometers  to  record  about  120  teleseismic  and  80  local  and  regional 
earthquakes.  The  broadband  seismometers  used  in  the  deployment  were  Guralp  CMG- 
3T’s.  The  data  were  continuously  recorded  by  Reftek  72-06  24  bit  digital  recorders  using 
50  samples  per  second. 

The  quality  factor  Qp  was  estimated  for  southern  Tibet  by  computing  spectral  ratios  of 
thirty  seconds  of  P-wave  train  in  narrow  frequency  bands  for  five  local  earthquakes.  Figure 
2  shows  seismograms  for  a  local  event  recorded  at  two  epicentral  distances  along  with  the 
correpsonding  spectra  for  the  P-wave  coda.  The  results  of  the  attenuation  calculation 
are  shown  in  Figure  3.  The  theoretical  curve  correspds  to  a  model  random  medium  with 
7  =  0.2,  background  velocity  cq  =  7.5  km/sec,  correlation  length  a  =  1.0  km,  and  fractal 
dimension  D  =  9/4.  Accordingly  the  high  frequency  dependence  of  Qp^  is 

~  (6) 

The  seemingly  large  value  of  7  could  be  reduced  by  including  intrinsic  attenuation  in  the 
calculation. 

Intermittent  Inhomogeneity 

Intermittency,  as  defined  by  Mandelbrot,  is  the  fact  that  in  natural  turbulence  dissi¬ 
pation  is  not  distributed  uniformly  in  space  (Mandelbrot,  1977).  Likewise,  there  is  no  a 
priori  reason  to  believe  that  seismic  velocity  heterogeneity  is  homogeneously  distributed 
in  three  dimensional  Euclidean  space.  Instead,  some  crustal  regions  may  be  very  hetero¬ 
geneous  while  others  are  devoid  of  velocity  impedence.  Furthermore,  the  region  where  the 
scalar  field  quantity  concentrates  may  sufficiently  tortured  and  convoluted  such  that  the 
geometry  is  appropriately  described  by  a  fractal.  The  case  of  homogeneous  turbulence 
concentrated  in  a  space  of  dimension  jD  <  3  has  been  called  fractally  homogeneous  turbu¬ 
lence  by  Mandelbrot  (Mandelbrot,  1976).  Likewise  the  case  of  seismic  velocity  fluctuations 
concentrated  on  a  fractal  may  be  called  fractally  homogeneous  inhomogeneity.  Thus  there 
are  two  dimensions  which  describe  general  fractally  homogeneous  fields.  The  dimension  D 
describes  the  way  the  field  is  spread  around  while  the  dimension  D  describes  the  way  the 
field  is  connected.  An  important  result  of  intermittency  is  that  the  spectra  of  fluctuations 
is  altered  by  a  factor  related  to  the  dimension  of  the  support  of  the  fluctuations.  The  sup¬ 
port  of  a  function  is  the  set  of  values  of  the  axgument  for  which  the  function  is  non-zero. 
For  example,  in  the  case  of  fluid  turbulence  where  the  structure  function  goes  like  the 
spectrum  is  proportional  to  where  B  is  related  to  a  function  that  describes  the 

process  of  ‘curdling’  which  concentrates  the  inhomogeneity  (Mandelbrot, 1976).  Likewise 
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for  fields  with  structure  functions  like  the  spectral  exponent  is  —(2(3  —  D)  +  B  +  l) 

so  that  within  an  intermittent  version  of  the  model  described  above  the  scattering  atten¬ 
uation  quality  factor, 

~  k-^2(3-D)+B)  (7) 

In  the  case  of  absolute  curdling,  the  support  of  the  field  is  a  closed  set  and  the  parameter 

B=-{3-D){2{Z-D)-l).  (8) 

If  the  dimension  of  the  support  D  lies  between  2  <  D  <  3  then  0  <  B  <  2/3(2(3  —  D)  —  1). 
As  an  example,  consider  a  random  fractal  medium  with  D  =  2.25  which  has  curdled  into 
a  support  with  fractal  dimension  D  =  2.5.  One  finds  that  B  =  —If A  and 

~  (9) 

In  general,  if  the  model  velocity  field  is  indeed  intermittent  {D  <  3)  and  and  the  iso¬ 
surfaces  have  dimension  D  >  2.5,  the  factor  B  is  positive.  This  implies  that  if  the  observed 
frequency  dependence  of  scattering  Q  is  attributed  totally  to  homogeneous  inhomogeneity, 
the  amount  by  which  the  iso-surface  dimension  exceeds  two  would  be  overestimated.  On 
the  other  hand,  if  D  <  2.5  the  factor  B  is  negative  and  homogeneous  theory  yields  an 
underestimate  of  the  dimension  D  of  constant  velocity  surfaces. 

Several  studies  of  coda  wave  attenuation  have  yielded  frequency  dependent  quality 
factors  Qc-  Some  results  are  summarized  below  cis  a  function  of  tectonic  region.  The  value 
of  the  iso-surface  dimension  D  is  given  for  three  different  dimensions  of  support  D. 


Region 

Western  U.S.  (Singh  and  Hermann,  1983) 
Garwhal  Himalaya  (Gupta  et.  ai,  1995) 
Western  Anatolia  (Akinci  et.  a/.,  1994) 
Southern  Tibet  (Qp)  (This  study) 


D{D  =  3) 

D{D  =  2.75) 

D{D  = 

0.7 

2.65 

2.70 

2.80 

0.9 

2.55 

2.57 

2.60 

1.0 

2.5 

2.5 

2.5 

1.5 

2.25 

2.17 

2.0 

Interestingly,  when  /3  =  I  intermittency  has  no  effect  on  the  fractal  dimension  D. 


Conclusions  and  Recommendations 

Scattering  attenuation  measurements  provide  an  effective  means  to  constrain  the  an¬ 
alytical  and  geometrical  properties  of  seismic  media  with  random  velocity  perturbations. 
Within  a  fractal  random  media  model,  the  attenuation  calculation  for  southern  Tibet  yields 
a  fractal  dimension  of  velocity  iso-surfaces  of  £>  =  9/4  which  implies  the  existence  of  com¬ 
plex  large-scale  structure  in  the  upper  lithosphere.  For  the  more  general  case  of  fractally 
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homogeneous  heterogeneity,  the  frequency  dependence  of  the  scattering  quality  factor  is  a 
function  of  two  parameters  D  and  D.  The  predicted  frequency  exponent  can  differ  from 
the  homogeneous  case  by  as  much  as  2/3  for  reasonable  choices  of  the  parameters  D  and 
D.  However,  the  relative  contributions  of  topological  disorder  and  intermittency  are  inex¬ 
tricably  linked  and  cannot  be  isolated  by  scattering  attenuation  measurements.  Continued 
investigation  of  wave  propagation  in  fractally  homogeneous  media  may  further  quantify 
the  statistical  properties  of  crustal  heterogeneity. 
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Seismograms  for  a  local  earthquake  recorded  at  two  epicentral  distances 
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Figure  2.  Vertical  traces  for  a  local  event  of  magnitude  4.7 
recorded  at  two  epicentral  distances.  The  corresponding 
P-wave  coda  spectra  are  plotted  on  the  bottom  demonstrating 
the  frequency  dependent  attenuation. 


480 


High  Frequency  P-wave  Quality  Factor 


Figure  3.  Observed  quality  factor  Q  for  P-wave  coda.  All 
the  data  for  five  local  earthquakes  in  Tibet  are  plotted  as  open 
circles.  The  mean  in  each  frequency  band  is  plotted  as  an 
asterisk  and  the  straight  line  is  fit  in  the  least  squares  sense. 
The  frequency  dependence  is  Q  "  f ’  ® 
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